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Abstract 

After summarizing briefly some numerical results for four- dimensional supersymmetric 
SU(2) Yang-Mills quantum mechanics, we review a recent study of systems with an infinite 
number of colours. We study in detail a particular supersymmetric matrix model which 
exhibits a phase transition, strong-weak duality, and a rich structure of supersymmetric 
vacua. In the planar and strong coupling limits, this field theoretical system is equivalent 
to a one- dimensional XXZ Heisenberg chain and, at the same time, to a gas of g-bosons. 
This not only reveals a hidden supersymmetry in these well-studied models; it also maps 
the intricate pattern of our supersymmetic vacua into that of the now-popular ground 
states of the XXZ chain. 
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1 Introduction 



This lecture reviews a recent progress in studying the large N limit of simple supersym- 
metric quantum mechanical systems which result from the dimensional reduction of field 
theories with gauge symmetry. Models of this type have been studied for a long time 
[HE], [3] and have many different applications, depending on the space-time dimension, D, 
of the unreduced theory [HE] (cf. Tabled]). For D=2 they can be often solved analyti- 
cally [H IU E] [7] , providing quantitative realization of supersymmetry. In four space-time 
dimensions one of them is nothing but the small volume limit of the Yang-Mills gluody- 
namics revealing the spectrum of zero volume glueballas as a special case |S1 [9] HOI CD] • 
Finally, for D=10 they make contact with the M-theory via the BFSS hypothesis |12j . 
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Table 1: A "road map" of a parameter space of supersymmetric Yang- Mills quantum 
mechanics. Black dots denote models well understood/solved, while open circles mark 
ones under study. A checkmark labels parameters of a system discussed in this lecture. 

Although we shall be mainly concerned with the large N limit, we would like to give 
one example of the N=2 model. 

1.1 Four-dimensional supersymmetric Yang-Mills quantum me- 
chanics at finite N 

The system has three bosons and two fermions, both in the adjoint representation of 
SU(2). The Hamiltonian reads [1] 

H = ^P l a p l a + ^e abc e ade xlx 3 c x l d xl + l -^-e abc iplT k ip b x k c , 
i = l,2,3 a =1,2,3. 

The spectrum is obtained numerically by diagonalizing H in the gauge invariant eigen- 
basis of the occupation numbers of all degrees of freedom n % a = o) a a l a ,rf a = p a f l a [3] 

\{<,4}>= e «w:/ t r/ t :...io>. 

contractions 
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The Hilbert space was cut by restricting the gauge invariant total number of bosonic 
quanta 

B = a\ a \ < B max . 

6,i 




Figure 1: The spectrum, and its supersymmetry structure, of the D=4 supersymmetric 
Yang-Mills quantum mechanics. 

It turned out that the physical (i.e. convergent with the cutoff eigenenergies) could 
be computed for sizes of bases well within a reach of a reasonably fast PC. The spectrum 
obtained in this way is shown in Fig. [T] [5J. It reveals a series of dynamical supermulti- 
plets with degenerate eigenenergies. The eigenstates within a supermultiplet are indeed 
the supersymmetric images of one another. These states are localized and they form a 
discrete spectrum. In addition there is also a continuum. It occupies central sectors of 
the figure, i.e. ones with the conserved fermion number F = 2,3,4. Supermultiplets of 
non-localized states occur at every second value of the angular momentum. Moreover, 
in these sectors the discrete and continuous spectra coexists at the same energies. This 
unusual feature of gauge interactions has never been observed so directly before. There 
are two supersymmetric vacua with fermion numbers F = 2 and 4. This is nothing but 
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the zero- volume manifestation of the existence of the AA condensate in the space extended 
theories [IB] . The condensate assumes two different (in fact opposite for SU(2)) values, 
again in agreement with the non-trivial predictions of the unreduced theory [TJ1 [15] . 

2 The large N limit and the planar calculus 

Above approach, even though quite successful, is naturally limited to the low number 
of colours. For higher N the Hilbert spaces become too large. Fortunately at infinite 
N Fock bases simplify enormously allowing again to reach quantitative results [T6]-|20j. 
This simplification is usually phrased in terms of the Feynman diagrams and topological 
expansion [21]. However it can be also formulated in the way suitable for the Hamiltonian 
formalism [TH1E]. To this end introduce matrix creation and annihilation (c/a) operators. 

a ik = V2a a T? k , f ik = V2f a T? k , i, k = 1, N. 

The physical basis can be conveniently generated by acting with the gauge invariant 
building blocks (bricks) 

(aW), (aWa^), (aWaW), . . . , (af N ), (.) = Tr[.) 

and their products, on the empty Fock state |0). Similarly for the fermionic sectors one 
employs bricks with mixed fermionic and bosonic operators. It turns out that in the 
large N limit only the single trace operators are relevant. All products of traces are 
either non-leading or do not provide new information (see [22] for recent results on that 
point). This observation culminates in the set of simple rules to calculate explicitly matrix 
elements of various Hamiltonians. One basically applies the Wick theorem employing the 
commutation rules 

[an,., a ji] = SuSkj 

and identifying the leading contributions. Details of such "planar calculus" have been 
presented in [161 El EEH] • Here we summarize only two examples 
A normalized state with n gluons in the F=0 sector reads 

|n >=^fr [(at)"] |0 > . (1) 

The normalization factor 

AC = < 0|Tr[a n ]Tr[(a t )"]|0 > 

= < 0|(12)(23)...(nl)[l / 2 , ][2 , 3'] . . . [n'l']|0 >, 
(12) = a hh , [12] = a ] hi2 . 

receives the maximal contribution only when the adjacent creation and annihilation op- 
erators are contracted. This gives 

AC = nN n . 
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Similarly a matrix element of a typical term in a generic Hamiltonian 

H n +2,n = g 2 < n + 2\Tr[a) a) a) a]\n > 
can be explicitly calculated 

H n +2,n = g 2 N^n(n + 2), 

and shown to depend only on the 't Hooft coupling A = g 2 N. 

The rest of this lecture will be devoted to a rather simple system which turns out to 
be surprisingly rich. The model may be considered as a distant cousin of the D — 1 + 1 
super symmetric Yang-Mills theory reduced to the one point in space. 

3 A simple supersymmetric system 

Consider one fermion and one boson with the following supersymmerty generators and 
the Hamiltonian |16j . 

Q = V2Tr\fa*(l + ga% Qt = ^2Tr[f{l + ga)a], 

H = {Q, Qt} =H B + H f . (2) 

Explicitly 

tt f / f 2\ 2 +^ 2 

Hb — a'a + g[a ] a + a'a ) + g a' a, 
Hf = ff + g{f ] f{a ] + a) + /t( a t + a)f) 

+ g 2 (pafJ + /taat/ + /t/ a t a + /tat/a). 

This Hamiltonian conserves the gauge invariant fermion number F = Tr[f'f] and can 
be diagonalized in sectors with well defined F. For any finite N calculation of matrix 
elements of H quickly becomes cumbersome. At infinite N, however, the planar rules 
illustrated above give simple and compact expressions. 

3.1 Hamiltonian matrix at large N 

The gauge invariant basis in the F = sector can be chosen as ([T|) with n = 0, 1, 2, 3, .... 
Only Hb contributes in this case. The matrix elements in the planar approximation read 

<0,n\H\0,n> = (1 + A(l - 5 n i))n, 

< 0,n+ l|F|0,n >= < 0,n\H\0,n+l >= y/\y/n{n + l). (3) 

In the sector with one fermion the basis is 

\n>=^ r Tr[(pa^r}\0>, n = 0,l,2,.... 
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Now, and for higher F, both Hb and Hp contribute. The Hamiltonian matrix is again 
simple 



The system reveals many interesting features already in these two lower sectors [IB] , 
therefore we postpone the discussion of higher F's and turn to the physics of ensembles 
with at most one fermion. 

3.2 Numerical results 

The Hamiltonian matrix is sparse but infinite with rows and columns labeled by the 
gauge invariant number of bosons B = Tr[a^a\. To obtain, the spectrum numerically we 
introduce the cutoff, B max , limiting the number of bosonic quanta 



and increase the cutoff until the spectrum converges. Results are shown in Figj2]for few 
values of the 't Hooft coupling. The convergence with B max is satisfactory for A 7^ 1 
and is faster for lower eigenvalues. There, the limiting, "infinite volume", results can be 
easily recovered. The less trivial way so see this is to check for the supersymmetry which 
is broken by the cutoff and should be recovered only at infinite B max . Figure [3] shows 
first few energy levels in four lower fermionic sectors. The degeneracies between bosonic 
and fermionic partners of the F = (0, 1) supermultiplets are excellent. Moreover, there is 
also an unbalanced supersymmetric vacuum state with zero eigenenergy 0. These results 
provide also the non-trivial test of our planar rules: the planar approximation does not 
break supersymmetry. 

The slowing of the cutoff dependence around A = 1 is a characteristic signature of a 
phase transition. At this point the system looses the energy gap and the spectrum becomes 
continuous. Any finite number of low-lying levels collapses to zero at infinite cutoff, but 
the cutoff dependence becomes characteristic of that for the continuous spectrum. 

Another interesting feature of this phase transition is the rearrangement of members of 
supermultiplets. This is shown in FigJH where the A dependence of the first four eigenen- 
ergies from both sectors is displayed for few values of B max . Away from the critical region 
supersymmetry is quickly restored and the F=0 and F=l levels are undistinguishable 
even for low cutoffs. Close to the criticality partners do not have the same energies since 
SUSY is broken at finite cutoff. Interestingly however, they rearrange while the 't Hooft 
coupling passes its critical value. In particular, the new vacuum state appears when A 
moves from the low to the large coupling phases. This also implies (and was readily found) 
that the Witten index (restricted to the F = 0,1 sectors) jumps by one unit across the 
phase transition point. 

1 An apparent lack of the degeneracy between some states with higher F will be explained later. 




(4) 



B < B, 
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1 = . 2 



1=1 . 




Figure 2: The cutoff dependence of the spectra of H in the F=0 sector, and for a range 
of A's 

Finally, the system has the exact strong-weak duality 

6 (Bf-"(l/f.) - i) = I (Er^i) - b 2 ) . 

It follows directly from the matrix representation (J3J) for F — 1 . However it is not obvious 
in the F = sector, cf. Eq. ([3]), but it is a direct consequence of supersymmetry. 
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Figure 3: Low lying bosonic and fermionic levels in the first four fermionic sectors. 



3.3 Analytic solution 

All above results have been subsequently derived analytically [16] . For example the second 
vacuum state has the form b = vA 



i<».=E(t)"> 

n=l v ' v 



n) 



(5) 



It is indeed annihilated by Q\ and exists only for b > 1. 

Surprisingly, one can find analytically the complete spectrum and construct all eigen- 
states in the F = sector. To this end it is convenient to introduce another, not or- 
thonormal basis, 



\B n ) = \fn\n) + bVn + l\n + 1) . 

This basis is convenient since the action of the Hamiltonian on \B n ) is so simple that 
the generating function for the expansion of the eigenstates in that basis 



fix) = J^c n x n 



J~]cn\B n ) 



n=0 



n=0 



can be constructed. The solution reads 

1 1 



fix) 



fix 



a x + 1/b 
1 1 



F(l, a: 1 + a: 



x + b 



1 — a x + b 
E = a{b 2 - 1), 



F(l,l-a;2-a; 



x + 1/b 

x+l/b 



b<l, 



x + b 



b>l, 



and together with the quantization condition 

/(0) = 0, 



(6) 
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Figure 4: Lowest bosonic and fermionic levels, as functions of A, for different cutoffs 



reproduces numerical eigenvalues obtained earlier in the infinite cutoff limit. 
As a one check, set a = in the generating function, for b > 1, to obtain 

which indeed reproduces the second ground state ((5]). 
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4 Higher fermionic sectors: F=2,3 

Single trace states with two fermions 

\n u n 2 ) = -?—Tr^ ni fa^f]\0), 

are labeled by two integers whose ordering is important modulo a cyclic permutation. 
Pauli principle eliminates states with n\ = n 2 , hence we can always take n 2 < n\. The 
planar rules give for the Hamiltonian matrix in this sector [19] 

(ni,n2|i?|ni,n 2 ) = (m + n 2 + 2)(1 + b 2 ) - b 2 {2 - 6 ni>0 ) - 2b 2 5 n2tni+1 , 
(ni + l,n 2 |-H"|rii,n 2 ) = b{n x + 2) = (n x , n 2 |^|ni + l,n 2 ), 
(ni,n 2 + l|if|ni,n 2 ) = 6(n 2 + 2) = (m, n 2 \H\rii, n 2 + 1). 
(ni + 1, n 2 - l|i?|ni, n 2 ) = (n 1 ,n 2 \H\ni + 1, n 2 - 1) 

= 26 2 (l-5„ 2 , ni+1 ). 

Planar states in the three fermion sector are labeled by three integers modulo cyclic 
translations 

\n h n 2 ,n 3 ) = ^—Tr[at ni fW»W*fi]\0), 
< ni, ni < n 2 , ni < ^3. 

Again the Hamiltonian matrix was explicitly calculated 

{ni,n 2 ,n 3 \H\ni,n 2 ,n 3 ) = (m + n 2 + n 3 + 3)(1 + b 2 ) 

—b 2 (3 — 5 nit o — 5 n2 fi — 5 n3t o), 

(ni + 1, n 2 , nsl-fflni, n 2 , n 3 ) = b{n 1 + 2)A = (ni, n 2 , n 3 \H\m + 1, n 2 , n 3 ), 
plus cyclic 

(ni +l,n 2 — l,n 3 \H\ni,n 2 ,n 3 ) = b 2 A = (rii,n 2) n 3 \H\ni + l,n 2 — l,n 3 ), 
plus cyclic. 

where A = l/y/3 if ni = n 2 = n 3 , and A = v^3 if the final state is of this form, otherwise 
A = 1. 

Numerical computation of the spectrum proceeds as before. There is again a phase 
transition causing the critical slowing down around A c = 1. Away from it, the eigenener- 
gies converge satisfactorily, with the infinite volume results shown in Figj3l The spectrum 
in higher fermionic sectors is essentially different from the F — 0, 1 case. We again find 
dynamical F = (2,3) supermultiplets, however not all states with three fermions have 
their superpartners in the F = 2 sector. Instead, they are degenerate with ones from the 
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F = 4 sector. This pattern continues now for all F as follows from counting the number 
of states: it is an increasing function of F. Another novel feature is that the spectrum 
is rather irregular, while the eigenenergies are almost equidistant for F = 0,1. There 
is again the rearrangement of members of supermultiplets across the phase transition. 
Interestingly however, the two new SUSY vacua with F = 2 appear in the strong coupling 
phase, while there is none at weak coupling. 

Similarly to the F = case, the two vacua can be constructed analytically. To this 
end consider the "extreme" strong coupling limit of the Hamiltonian (|2J)[19J 

H sc = lim \h = Tr(ff) + ±-[Tr(a^a 2 ) + Tr^faf) + Trtfctfa)]. (8) 

A^oo A iV 

Surprisingly, Hsc conserves the number of bosonic quanta as well, and proves very useful 
in mapping the structure of the model in a//fermionic sectors (cf. the following Section). 
Coming back to F — 2, one can identify the finite dimensional (F = 2, B) sectors with 
zero eigenvalues, and construct the corresponding eigenstates. They read 

\F = 2,v)T = |0,1), (9) 
\F = 2,v)? = |0,3)-2|1,2), (10) 

the superscript referring to the infinite value of the 't Hooft coupling. These states can 
then be used to construct the two vacua at finite A > 1 

oo 

|F = 2,^ = ^(-v / A)-"|0,n), 

71=1 

\F = 2, v) x 2 = (1 + Q S H~ X Q^ W )~ X \F = 2, v)?, (11) 

with Q s = gTr[fa)\ H s = {Q„Q\}, Q^ w = Tr[pa] [12]. 

The restoration of supersymmetry can also bee seen on the level of eigenstates. Mem- 
bers of supermultiplets transform among themselves by SUSY charges Q and Q\ which 
become the well defined matrices in the planar bases. To see SUSY in this way define the 
following supersymmetry fractions 

q mn = J 2 <F + l,E m \Qi\F,E n >, (12) 

V -frm + &n 

which are the coordinates of the supersymmetric images of eigenstates in the correspond- 
ing fermionic sectors. Figj5] shows that SUSY fractions indeed quickly stabilize with 
incrasing B max . 

Supersymmtery fractions are also useful in defining the restricted Witten index which 
smoothly interpolates, at finite cutoff, between the two phases. The straightforward 
restriction of the sum 
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Figure 5: First five supersymmetry fractions. 
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Figure 6: Behaviour of the restricted Witten index, at T = 6, around the phase transition. 



to the F = 2, 3 eigenstates does not suffice since some of the states with three fermions 
remain unbalanced. Instead, one can define 



W R (T, \)=Ys( e 



-TEi _ e -TEi 



J2f E f\qfi 



i.e. we take as the energy of the supersymmetric partner the energy weighted by supersym- 
metric fractions. This definition enforces summation only over complete supermultiplets 
away from the critical region, while provides the smooth smearing among possible candi- 
dates around A c . The index defined this way changes smoothly at finite cutoff, and varies 
by two units as expected (cf. FigE]) 
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5 Arbitrary F 



A simple generalization of the F = 2, 3 cases for arbitrary number of fermions reads [18] 
\n) = \m, n 2 ,..., n F ) = * Tr^/V^t . . . a^ F f)\0) (13) 

M{n} 

So the states of the basis may be labeled by F bosonic occupation numbers {rii, n 2 , . . .n F } 
(configurations) modulo cyclic shifts. Some states are excluded by the Pauli principle. 
For example 

{3,3}, or {2,1,1,1,2,1,1,1} 

are not allowed, since they change sign while retaining their identity after the suitable 
number of cyclic shifts. When calculating normalizations and matrix elements one should 
keep track of symmetry factors defined as the number of cyclic shifts which bring a state 
back to itself. For example the state |2, 1, 2, 1, 2, 1) has the symmetry factor d = 3. 
Alternatively one can label planar states with the periodic binary strings (necklaces) 
with, e.g., zeros/ones corresponding to the bosonic/fermionic creation operators (for a 
more complete discussion see [IB"]). 

Detailed calculation of spectra for arbitrary fermion number is now in progress. How- 
ever many properties, e.g. the supersymmetry structure, has been inferred from the strong 
coupling limit of the model. 

As already mentioned, the strong coupling Hamltonian (|S]) conserves both fermionic 
and bosonic numbers of quanta. Hence the Hilbert space splits into the finite dimensional 
sectors, cf. Table El and the Hamiltonian in each sector becomes a finite matrix. Still the 
Hsc has the exact supersymmetry with the strong-coupling charges 

Qsc = -^ Tr if aj \ Q'sc = -^ Tr [/ flt2 ]' ( 14 ) 



acting along the "diagonals" 2F + B = const, of Table [2j Diagonalizing few of these 
matrices we have indeed found expected degeneracies^]. Interestingly, we have also found 
experimentally [TBI [20] that the zero energy eigenstates are located only in the sectors 
with 

B = F±1, F even. (15) 

which form the regular "magic" staircase, shown in the boldface, in Table [2j 

This observation, when combined with the weak coupling (harmonic oscillator) limit, 

explains the structure of SUSY vacua for all F and at any value of the 't Hooft coupling. 

Namely, for any even F there are two supersymmetric vacua for any A in the strong 

coupling phase. On the other hand, there is only one SUSY vacuum in the weak coupling 

phase and it has F = 0. 

The question why the supersymmetric vacua are located in the magic sectors (IT51) will 

be answered in the next Section. 



2 Obviously no cutoff was required this time. 
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Table 2: Sizes of gauge invariant bases in the (F, B) sectors. 

6 Two equivalencies with statistical systems 

Interestingly, our strong-coupling model is equivalent to the two well known and nontrivial 
statistical systems [20J. 

6.1 A gas of q-bosons 

Consider a one dimensional, periodic lattice with size F. At each lattice site put a bosonic 
degree of freedom described by its c/a operators a^, aj. The strong coupling Hamiltonian 
dU) is equivalent to the following one expressed solely in terms of bosonic variables 

F F 

H = B + Y, <W + + fri&Vi, (16) 

1=1 i=l 

where Ni = a\a>i, B = N\ + N2 + . . . + Np, and the action of bi, b\ is 

6t| n ) = | n + i) > b\n) = \n-l), 6|0) = 0, (17) 

that is, they create and annihilate one quantum without the usual y/n factors. They 
satisfy the following commutation rules 

[6,61=^,0 (18) 

with the same non-linear (in b's) 5 operator as in the Hamiltonian. Using the planar rules 
discussed earlier, one can show that the action of the first two terms of (jHJ) on the planar 
basis is exactly the same as the action of the first two terms of (TIB"]) in the eigenbasis of 
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Ni [20] • The last two terms of (|S]) are the same as the hopping terms of f)16p . Since the 
planar states acquire a phase (— upon cyclic shifts, the planar system has the same 
spectrum as (flfij) in the eigen-sector of the lattice shifts, U, with = (— l)^" 1 . This 
equivalence was cross-checked numerically for 3 < F, B < 7. 

The c/a operators are known in the literature [231 121j- Their algebra is the special 
case of the q-deformed harmonic oscillator algebra 

[b,tf}=q- 2N (19) 

with q —>■ oo. Transitions ffTTj) without the square roots are often referred to as assisted 
(i.e. independent of the occupancy) transitions. The system of q-bosons described by the 
Hamiltonian ([TBI is, to our knowledge, regarded as non-soluble. In view of the second 
equivalence, discussed below, it turns out to be in fact soluble. 



6.2 The XXZ Heisegberg chain 

The second equivalence has been proved employing another representation of the planar 
states ([I3D [20] 

-L\Tria^ mi ^r^ m2 (fr ■ ■ ■ (PT k )\0) = |(0) mi (l) ni • • • (0)^(1)"-} r > 1, m^, > ($0) 

Where the states are labeled by the binary strings with, e.g. 0/1 corresponding to the 
bosonic/fermionic creation operators. In this basis, the action of the strong coupling 
Hamiltonian ([8]) is equivalent to that of the XXZ chain 



i=l 



for the particular values of the anisotropy parameter A. The detailed correspondence 
reads 



r(+V2) 
XXZ 

-Hl~M 2) + \L , F even, B odd . 



H sc (F,B) = {- H m 2) + f T > n F ° d ^ (21) 



where the spin Hamiltonians are restricted to the translationally invariant sectors with 
the lattice size and the total spin fixed by (F, B): 

L 1 

L = F + B, S* = J3a? = -(F-B), (22) 
i=i 

Remarkably, this equivalence implies the existence of the magic staircase of the super- 
symmetric vacua discussed above. More than thirty years ago Baxter has found that, 
for A = —1/2, the ground states with S z = ±1/2 have particularly simple eigenenergy 
Eq = — \L for infinite L (25]. Recently his findings have been extended by Riazumov and 
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Stroganov to any finite, odd L [2B]. In view of (j21||22|) the Riazumov-Stroganov states are 
nothing but our strong coupling vacua [20J. Moreover, supersymmetry requires a host of 
degeneracies among the, seemingly unrelated, excitations of the XXZ chain. Amusingly, 
since the supersymmetry transformations change the lattice size, these excitations live on 
different lattices! [f| 

Finally, since the XXZ model is exactly soluble (e.g. by the Bethe Ansatz [28J), our 
super symmetric system is also soluble at infinite coupling. As one application we give 
here the algebraic determination of the Bethe phase factors for the first three steps of the 
magic staircase. The eigenenergies of Hxxz(A) are [29] 

A m 
E XXZ (A) = -L- + 2mA - 2 ^ cos Pi , (23) 

where the momenta —n < pj < it satisfy the following set of Bethe equations 

J?L. pi-Pl _|_ p-i-Pj _ OA 

e lL ^ = (-1)"'- 1 FT e'fo-pO e ^ e . tf± j = 1, ... ,m. (24) 

v ' J- J- pWj + e -m — 2 A 

1=1 

With m denoting the number of down spins in a chain. For the supersymmetric model 
m = B and (j23J translates into 

E SC (F } B) = F + 2 J2f = i cos Pjl for F odd , A = +i, (25) 

£ 5C (F, B) = F — 2 J2f =1 cospj, for F even, and 5 odd , A = -J. (26) 

J 2 

Solving numerically Eqs. (l24|) we have found that the supersymmetric vacua occur at F 
even and B odd, and are always given by the yet simpler sub-ansatz of the Bethe Ansatz 

Pi = 0, p 2 k+i = -P2k, fc = l,...,(B-l)/2. (27) 

which reduces the number of independent variables roughly by a factor of 2. Still, the 
Bethe equations can only be solved numerically in general. However, for the first three 
sectors the problem can be managed algebraically. We shall discuss it separately for each 
sector. 

6.2.1 F=4, B=3 

This corresponds to L = 7, m = 3, and Eq. (1271) implies 

Pi = 0, p 2 = p, p 3 = -p, (28) 

3 Supersymmetry of the, A = —1/2, XXZ chain and other statistical systems, has been discussed in 
the literature, albeit in slightly different contexts (see, e.g. [27] and references therein). We thank Jan 
de Gier for bringing this to our attention. 
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and f[2"4"j) reduces to (z = e ip ) 

z 6 = 1. (29) 
The admissible solution, which gives Esc — 0, is 

z = e'*. (30) 

6.2.2 F=4, B=5 

Now L = 9, m = 5, and with the aid of (|2"7j) Bethe equations reduce to (x = e ip2 , y = e WA ) 

xy + x + y xy + x + 1 



x + y + l xy + y + 1' 



y8 = x ^ + x + V xy + y + 1 ^ 
x + y + lxy + x + 1 

Solving (I3l~j) is difficult in general, however one can easily find solutions with vanishing 
energy (|2"6"j) . 

E 5C (4,5) = 2-2fa; + y + - + -) =0. (32) 
V x yj 

Introducing two symmetric variables 

s = x + y, p = xy, (33) 
one obtains from fl3T|) . (1321 the following two equations for s and p 

±p\ (34) 



8 + 1 

s + - = l. (35) 
P 



The admissible solution is on the negative branch and reads 



x 

64 



l - ^16 + iV2\J 15 + V33(7 - v 7 ^) - 4^-16(3 + v 7 ^) - ^v 7 ^ 15 + v^3(9 + V33) 
^ ^16 + iV2\/ 15 + VS3(7 - V33) + 4^-16(3 + V33) - i2^2\J 15 + ^(9 + V33) 



This pair indeed satisfies both Bethe equations, together with (I3"2"j) . and therefore corre- 
sponds to our supersymmetric vacuum. 
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6.2.3 F=6, B=5 

This sector corresponds to L = 11, m — 5. Now the reduced Bethe equations are 

10 xy + x + y xy + x + 1 

x = , 

x + y + 1 xy + y + 1 

10 xy + x + y xy + y + 1 

y = — i — n , — rr- 37 

As previously, we look for the solutions which satisfy Esc(6, 5) = 0, that is 

11 s 2 

x + y + - + - = 2, or - = - . 38 

x y p 1 + p 

The second equation follows from the product of the Bethe equations 

P = ±7— > or - = -——jr. (39) 
1 + s p 1 =F P 

These equations reduce to the fourth order polynomial equation for p. Again there is only 
one admissible solution and it lies on the negative branch 



X 72 



l - ^36 + 11 + Vl3(7 + VlS) - Q^2\J 6(-3 + Vl3) + 11 + Vl3(-5 + j , 
^36 + iV2y/ll + Vl3(7 + v 7 !^) + 6(-3 + Vl3) + iyftyjll + Vl3(-5 + Vl3) j . 



(40) 



As the last exercise let us show that these numbers are indeed unimodular which is 
not obvious at first sight. However it is easy to calculate 

\x + y\ 2 = ^(7 + Vl3), \x-y\ 2 = ~(5 - Vl3), therefore \x\ 2 + \y\ 2 = 2. (41) 

Further, since x = a + b and y = a — b, one can observe that 

(ab*) 2 = ^(-11 + V13), => \x\ 2 = \y\ 2 , (42) 

which together with f j4TT) implies that x and y are indeed pure phase factors. As a 
byproduct one obtains 

Re{xy*) = cos (p 2 - p 4 ) = 1(1 + >/l3). (43) 

o 
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7 Discussion 



The direct diagonalization of a Hamiltonian matrix is usually considered a viable tool, 
for finding a spectrum, only for finite matrices. It turns out however, that the approach 
works in many cases with the infinite dimensional Hilbert space as well. Although such 
Hilbert spaces appear already in many quantum mechanical systems with finite number of 
degrees of freedom, the true challenge for the Hamiltonian formalizm is posed by the field 
theoretical problems. An interesting family of systems results from the dimensional reduc- 
tion of field theoretical models to a one point in space. As such, they again have the finite 
number of degrees of freedom, however they inherit many advanced features of parent field 
theories, for example their symmetries including supersymmetry. The straightforward di- 
agonalization proved quite successful in uncovering quite rich and nontrivial spectra of 
super symmetric Yang-Mills quantum mechanics in various dimensions. 

This lecture reviews the recent study of a system with an infinite number of degrees of 
freedom, namely a particular supersymmetric gauge quantum mechanics (also referred to 
as a matrix model) with the infinite number of colours. The model was conceived as the 
illustration of the planar calculus in the Fock space. However it was subsequently found 
that the system has an interesting physics which connects to many recently discussed 
issues. For example, the system undergoes the discontinuous phase transition in the 
't Hooft coupling which is accompanied by the remarkable rearrangement of dynamical 
supermultiplets. It enjoys the strong-weak duality in the lowest (and simplest) fermionic 
sectors where the complete, exact spectrum can be found analytically. The full structure 
of intervening supermultiplets begins with two fermions and goes on ad infinitum. In each 
bosonic sector two new supersymmetric vacua appear in the strong coupling phase while 
there is only one in the weak coupling region. The behaviour in the strong coupling phase 
has been found by studying the system at the infinite value of the 't Hooft coupling and 
then extending it to the whole strong coupling regime. 

At infinite coupling the model reveals also its connections with the statistical physics, 
which proves, among other things, that a quantum mechanics at infinite N becomes a 
bona fide field theory. The supersymmetric planar model is exactly equivalent to the 
one-dimensional quantum XXZ Heisenberg chain, and at the same time, to the one- 
dimensional lattice gas of q-bosons. The strong coupling vacua found in the SUSY matrix 
model turn out to be nothing but simple ground states of the XXZ chain which were 
found more than thirty years ago. In addition, the XXZ chain appears to have the hidden 
supersymmetry which results in many degeneracies among various energy eigenstates. 
An unusual feature of these supersymmetry transformations is that they connect states 
propagating on different lattices. 

Second equivalence is with a gas of q-bosons, with the infinite deformation parameter. 
The latter mapping holds exactly in all fermionic and bosonic sectors while the former 
does not work for even F and B. 

The chain of correspondences discovered in [20J implies also that the specific, nonlinear 
oo-bosonic Hamiltonian is in fact soluble via solubility of the XXZ chain. Vice versa: the 
same solubility implies that the supersymmetric matrix model is exactly soluble at the 
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infinite value of the 't Hooft coupling. It is also quite conceivable that this solubility holds 
in the whole strong coupling phase [T9~] . 

Summarizing: planar calculus applied directly in the Fock space turned out be rather 
promising tool in studying some simple, but non-trivial, models. It allowed to find phe- 
nomena which bear a tantalizing similarity to ones found in the much more advanced 
systems [30j-[35j. It remains to be seen if this approach can be extended to the full, i.e. 
space extended field theories. 
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